Absorbing Phase Transitions of Branching- Annihilating Random Walks 
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The phase transitions to absorbing states of the branching-annihilating reaction-diffusion processes 
mA — > (m + k)A, nA — * (n — I) A are studied systematically in one space dimension within a new 
family of models. Four universality classes of non-trivial critical behavior are found. This provides, 
in particular, the first evidence of universal scaling laws for pair and triplet processes. 
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Even though it is widely believed that far from equilib- 
rium phase transitions can be classified, recent numerical 
results reveal that we are still far from a satisfactory un- 
derstanding of even the relevant ingredients deciding to 
which class a given transition belongs. This is in partic- 
ular true for transitions from a fluctuating phase to one 
or several absorbing states (APT, for absorbing phase 
transitions) where, in spite of a wealth of analytical and 
numerical studies, rather little is known beyond the ex- 
istence of the prominent universality class of directed 
percolation (DP) [Jl|. In terms of reaction-diffusion pro- 
cesses, the DP class is often represented by the simplest 
branching-annihilating reactions: A —> 2A, A — > 0. 

Following an early suggestion by Grassberger a 
lot of attention has been devoted recently to the case 
of binary or pair reaction-diffusion processes such as 
2 A — > (2 + k)A, 2 A — > or A, which involve two parti- 
cles both for branching and annihilating || |], 0, ^. [/], 0, 
|, [H], [H| [0| [TJ, [lf§. Usually designated under the 
name of "pair contact process with diffusion" (PCPD), 
it is still unclear as of now if systems of this type can 
exhibit universal scaling laws or even if (simple) scaling 
occurs at all at the APT. The origin of the strong devia- 
tions to scaling observed is debated , and various 
conclusions have been drawn ranging from diffusion-rate 
dependent sub-classes (6| to continuously-varying expo- 
nents p3 to no scaling to slow crossover to the DP class 
Jl5| |. At the analytical level, a bosonic field theory of 
PCPD processes exists but is not renormalizable whereas 
no fermionic version is available Completing 
this unsatisfactory picture, similar results were reported 
recently for the "triplet contact process with diffusion" 
(TCPD) 3A -> 4A, 3A -> 2A, where three particles are 
involved for both reactions |l6| . Finally, the role of the 
conservation of parity of the number of particles in APT 
is not well understood either: it is known to be relevant 
for some one-dimensional DP-like processes, giving rise 
to the so-called "parity-conserving " (PC) class 0, ||, 
whereas it was argued not to influence the critical prop- 



which we draw a considerably clarified picture of the 
above situation. This is achieved thanks to a new class 
of very simple models which, contrary to all the works 
cited above, abandon the "fermionic constraint" usually 
considered for both theoretical and practical reasons (to 
avoid the divergence of the density of particles in the 
active phase, and — hopefully — increase numerical effi- 
ciency). Our results include the first evidence of univer- 
sal (ordinary) scaling for PCPD and TCPD processes, 
the classification of "hybrid" rules (i.e. those for which 
to 7^ n), and further insights into the conditions under 
which the conservation of the number of particles modulo 
2 or 3 is able to change the "reference" scaling laws. 

Following Hinrichsen JT(| , we may write, for the mA — > 
(m + k)A, nA — > (n — I) A processes, the following "mean- 
field Langevin equation" expected to govern the coarse- 
grained local density p: 



d t p = ap m - bp n - cp m+1 + DV 2 p + C(x, t) 



(1) 



erties of PCPD systems |19 

In this Letter, we report on a systematic numerical 
investigation of the reaction-diffusion processes mA — » 
(m + k)A, nA — > (n — I) A in one space dimension, from 



where a, b, and c are positive constants related to the 
reaction rates and D is the diffusion constant. Note that 
of the two negative terms in Eq. (|l]) only the lowest-order 
one is relevant. As usual, the stochastic nature of the 
process is embodied in £, a delta-correlated noise whose 
amplitude is a power of p (thus ensuring the absence of 
fluctuations in the p = absorbing phase) : 

(C(x, *)C(x', f')) = T/Ax, t)5 d (x - x!)6(t - f) . (2) 

In the absence of branching (i.e. in the inactive phase 
leading to an absorbing state), the noise dimension p = n 
and the upper critical dimension of the annihilation pro- 
cess is d c = 2/(n — 1). In the critical region, fluctuations 
are expected to increase the noise strength and thus to 
reduce p and increase d c . In the following, we restrict 
ourselves to m, n < 4 (in one dimension, non-trivial APT 
are mostly expected to occur if n < 4). 

Each of the branching-annihilating processes mA — > 
(to + k)A, nA — > (n — l)A is defined by the 4 integers 
(m,n,k,l). Obviously, one must have < n — I < m 
to insure the existence of at least one absorbing state. 
For legibility, m and n will be coded below by the letters 
s, p, t, q (for singleton, pair, triplet, quadruplet). The 
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FIG. 1: Time decay of the order parameters at criticality 
for the PCPD rule ppl2 (p c = 0.795410(5]) in logarithmic 
scales. Single run on a system of 2 22 sites, starting from 2 
particles on each site (other initial conditions do not influence 
the asymptotic behavior). Different order parameters, from 
top to bottom: density of particles, density of pairs (e.g., by 
summing all particles of sites occupied by at least 2) , fraction 
of occupied sites, and fraction of sites occupied by at least a 
pair. Inset: same, but for order parameters multiplied by t 02 . 



PCPD rule 2A -> 3A, 2A -> is thus noted ppl2, and 
ttxx denote TCPD processes. The fermionic constraint 
adopted in most PCPD and TCPD models studied so 
far can be seen as counter-productive: the actual imple- 
mentation of their Monte-Carlo simulations often lead to 
complicated rules, inefficient for both code elaboration 
and simulation. More importantly, the fermionic con- 
straint may well be at the origin of the strong deviations 
to scaling observed. Our branching-annihilating random 
walk (BARW) models are designed to bypass both of 
these problems. Particles of a single species A evolve in 
parallel in two synchronized sub-steps: random walk on 
the lattice (diffusion) followed by on-site reaction. For 
simplicity, in the following, the diffusion constant is kept 
constant: all particles always jump to a randomly chosen 
nearest-neighbor. Let us consider, to describe our on-site 
reaction scheme, the PCPD rule ppl2. (The generaliza- 
tion to all other rules studied here is straightforward.) 
Suppose that Ua particles are present at a given site. If 
Ua = 1 ; nothing happens (this is the main PCPD con- 
straint). For ua > 2, each of the \jia/2\ pairs into which 
the local population can be divided branches with prob- 
ability pL«A/2j (thug creating each time one new parti- 
cle for this particular example) otherwise it annihilates. 
The only parameter is p: for large p, branching is likely, 
and one expects to be in the active phase. For small p, 
annihilation dominates (indeed this is the only process 
at play for p = 0), leading to an absorbing state. One 
key feature is the nonlincarity introduced by raising p 
to the power \jia/2\: branching/annihilation is inhib- 
ited/enhanced for large local populations, preventing the 
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FIG. 2: Time decay of the total density of particles at crit- 
icality for various TCPD rules in logarithmic scales. Single 
runs on a system of 2 22 sites, starting from 3 particles on 
each site. From top to bottom: rule ttll, p c = 0.53948(1); 
rule ttl2, p c = 0.72899(1); rule ttl3, p c = 0.822055(3); rule 
tt22, p c = 0.61900(2); rule tt33, p c = 0.70323(2); Inset: same, 
but multiplied by t 0,27 . All data have been shifted for clarity. 



divergence of population of usual bosonic models. All 
the results presented below were obtained for reaction 
schemes of this type, but we have checked that the func- 
tional form of the nonlinearity as well as other details do 
not change the critical properties observed. 

As recommended when studying APT numerically, we 
first investigate, for a given rule, the decay of the or- 
der parameter from some highly active, correlation-less, 
initial condition in a large system. Above threshold, ac- 
tivity eventually reaches a constant level in time. Well 
below threshold, one observes a decay typical of the an- 
nihilation process (exponential for n = 1, algebraic for 
n > 1). At p = pc, in the usual framework, one expects 
a non-trivial algebraic decay characterized by the scal- 
ing exponent <5 = @/v\\- Once the threshold determined, 
other scaling exponents can be estimated. In this work, 
we also give — without showing the data — our estimates 
of z, obtained by the finite-size scaling of the mean life- 
time of the system at threshold ((t) ~ L z ), and of j3, 
given by the decay of the stationary density with the dis- 
tance to threshold (lim^oo p ~ (p — p c ) )■ More precise 
estimates will appear elsewhere ^0|. In the following, 
the role of conservation (modulo 2 or 3) of the number 
of particles is discussed at the end. Before proceeding to 
higher-order rules, we report that ssxx rules are easily 
verified to exhibit DP-class scaling laws (not shown). 

In Fig. [l] we show the critical behavior of the PCPD 
rule ppl2. Quite easily, a clean algebraic decay in time 
of all order parameters is observed over many decades 
with the exponent <5 PCPD = 0.200(5). At the determined 
threshold, the mean lifetime scales nicely with system 
size, yielding the estimate z PCPD = 1.70(4). Finally, the 
decay of the stationary order parameter with distance to 
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FIG. 3: Same as Fig. y but for hybrid rules m > n. From top 
to bottom: rule tpl2, p c = 0.81156(1), 2 24 sites; rule ts21, 
p c = 0.62 1 21(2), 2 20 sites; rule tp22, p c = 0.756625(5), 2 22 
sites; rule psll, p c = 0.71529(1), 2 20 sites. Inset: same, but 
multiplied by t 016 . 
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FIG. 4: Same as Fig. |2| but for ptxx hybrid rules. From 
top to bottom: rule pt42, p c = 0.56775(4); rule pt22, 
p c = 0.54640(2); rule ptl2, p c = 0.33250(3); rule ptl3, 
p c — 0.66210(2). Inset: same, but multiplied by t ' 2 . 



TABLE I: 


Critical exponents 


for basic classes 




Class 


S 


z 





DP (from Q) 


0.1595 


1.58 


0.2765 


PCPD 


0.200(5) 


1.70(5) 


0.37(2) 


TCPD 


0.27(1) 


1.8(1) 


0.90(5) 


PC (from H) 


0.286 


1.76 


0.92 



threshold allows to estimate /3 PCPD = 0.375(10). These 
exponent values appear to be universal: this has been 
checked to the above numerical accuracy for rule ppll 
and with less care for several other rules. 

Similar results were obtained for TCPD rules. Al- 
though scaling usually sets in later than for the PCPD 
rules, it is well established over a large range of scales. 
Again, scaling laws are found to be universal within nu- 
merical accuracy. Fig. |^ shows typical results obtained 
for the decay of the particle density at threshold for five 
different rules, from which we estimate £ TCPD — 0.27(1). 
Other exponents are reported in Table B. 

Next, we consider hybrid rules for which to > n. At 
the mean-field level, they are expected to present first- 
order transitions (see Eq. ([!])). Nevertheless, in one space 
dimension, all the cases we considered show DP scaling, 
at least for the decay exponent 8. This is true for n = 1 
(simple radioactive decay), a case where we have stud- 
ied rules psll, ts21, and qsll. But this is also true for 
higher-order annihilation processes, as shown by pair- 
annihilation rules tpl2 (2 variants) and qpl2, as well as 
by triplet-annihilation rule qtl2. As seen in Fig. ||, for 
high-order branching process scaling sets in rather late, 
probably because the initial conditions chosen are not 
"optimal" . That a given process for which a first-order 
APT is predicted at mean-field level exhibits a continuous 



transition in low dimensions is by no means surprising, 
especially in one dimension, where this has been noticed 
early pi| . Nevertheless, to our knowledge, satisfactory 
analytical arguments for our finding are not available. 

Hybrid rules for which the branching process is of lower 
order than the annihilation reaction (m < n) are also 
easily investigated within our family of models. First, we 
note that whenever n = 4, no non-trivial APT was found 
in any of the rules sqxx, pqxx, tqxx or qqa;x that we have 
considered. In other words, p c = and the decay is then 
that of the mean- field prediction (i -1 / 3 ), in agreement 
with the fact that the upper critical dimension of the 
decay process is | < 1. For to < n < 4, we find that the 
universality class of the critical behavior is determined 
by the order of the branching process: m = 1 rules such 
as spl2, stl3, st23 exhibit DP-scaling, in agreement with 
the analytical arguments of ]H| ]. Pair-branching rules 
(to = 2) ptl2 and ptl3, after a rather long crossover time, 
show scaling compatible with the PCPD class (Fig. ||). 
That to sets the universality class in this case is actually 
not too surprising: it is clear at the mean-field level of 
Eq. (Q) , and can also be deduced from "vertex generation 
arguments" such as those developed in p8[ . 

We now pay particular attention to the rules for which 
the total number of particles is conserved modulo 2 (par- 
ity conservation) or modulo 3. Phase space is then di- 
vided into 2 or 3 disjoint sectors, each of which may or 
may not possess an absorbing state. For example, the 
even sector of parity-conserving rule sp22, the archetyp- 
ical rule of the PC class, includes an absorbing state 
(the empty configuration) whereas its odd sector does 
not. For parity-conserving PCPD rule pp22, on the other 
hand, each sector has an absorbing state. We studied at 
least one rule in each of the relevant sub- families (namely 
pp, tt, sp, pt, tp, qp, qt for parity conservation, and st, 
pt, tt, qt for conservation mod 3). We find that when- 
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TABLE II: Universality class of BARW processes. "0" stands 
for no non-trivial APT (p c — 0). The second markings indi- 
cate the change of class whenever mod2 or mod3 conservation 
plays a role. See text for details. 



m\n 


1 


2 


3 


4 


1 


DP 


DP/PC 


DP/0 





2 


DP 


PCPD 


PCPD/0 





3 


DP 


DP 


TCPD 





4 


DP 


DP 


DP 






ever every sector includes an absorbing state, mod2 or 
mod3 conservation does not change the "reference" scal- 
ing class. Thus rule pp22 is in the PCPD class (as sug- 
gested in rules tt22 and tt33 in the TCPD class 
(Fig. g). Interestingly, rules tp22, qp22, qt22, and qt33 
are found to be in the DP class, like their non-conserving 
counterparts, and thus constitute the first examples of 
DP-class rules with mod2 or mod3 conservation (Fig. ||). 

Conversely, conserving rules with at least one non- 
absorbing sector do not offer such a clear conclusion: as 
expected, rules sp22 and sp42 fall into the PC class (not 
shown) . As argued in , rule st33 does not have a non- 
trivial APT. This is also observed for rule pt33, prob- 
ably for similar theoretical reasons. The only possible 
case left is that of parity-conserving hybrid rules pt(2/c)2. 
There the even sector is non-absorbing because the empty 
state — which is absorbing — cannot be reached from any 
other configuration. Numerical simulations of rules pt22 
and pt42, however, do not show any significant difference 
from their non- conserving cousins (Fig. ^). This result, 
which needs confirmation due to the lateness of scaling, 
seems to refute the suggestion of jl9| that influence of 
parity conservation may be equivalent to having a non- 
absorbing sector of phase space. Instead we would like 
to propose that rules sp(2/c)2 are the only ones in the PC 



that the critical behavior of all processes considered here 
is as follows (Table ||) : rules with four-particle annihila- 
tion (n — 4) do not have non-trivial APT, indicating (in 
partial agreement with [jl6|) that d c < 1 in this case. Hy- 
brid processes fall into one of these basic classes: m> n 
rules exhibit DP critical behavior, whereas the class of 
to < n rules is set by m. Finally, we have shown that 
mod2 (parity) or mod3 conservation seems to act only 
on sp(2/c)2 rules, suggesting that the PC class should be 
seen as the (one-dimensional, generalized) Voter class. 

Obviously, our numerical findings need to be confirmed 
by analytical approaches. We hope that, at the very 
least, they will trigger some new lines of attack to the 
difficult issues at play. At the numerical level, ongoing 
work aims at obtaining more comprehensive results, in- 
cluding the study of spreading exponents and the case of 
higher space dimensions. 

We thank F. van Wijland for useful discussions. 



class because they are the only BARW processes equiva- 
lcnt to generalized twb-statc Voter models, as defined in 
pjj , i.e. spin models with up/down symmetry, no bulk 
noise, and an order/disorder transition. This remark, 
which will be developed in p5[ , calls for renaming the 
PC class the Voter class. 

We now summarize and discuss our results. We first 
note that our BARW models, at criticality, exhibit or- 
dinary scaling after some crossover scale, and do not 
seem to be plagued by the strong deviations observed 
with even the simplest fermionic PCPD or TCPD rules 
p3| , p"5| . Numerical results of basic fermionic models can 
be shown to converge to the critical behavior found here 
ppf . To our numerical accuracy, the scaling exponents 
recorded at criticality lead to conclude to the existence of 
three basic universality classes: DP, PCPD, and TCPD. 
Table | summarizes our current estimates of the basic 
exponents 5, z, and /3, pending more precise ones pOj . 
Extrapolating from our numerical findings, we believe 
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